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ABSTRACT 

We numerically investigated the dynamical architecture of 47 UMa with the 
planetary configuration of the best-fit orbital solutions by Fischer et al. We 
systematically studied the existence of Earth-like planets in the region 0.05 AU 
< a < 2.0 AU for 47 UMa with numerical simulations, and we also explored the 
packed planetary geometry and Trojan planets in the system. In the simulations, 
we found that "hot Earths" at 0.05 AU < a < 0.4 AU can dynamically survive at 
least for 1 Myr. The Earth-like planets can eventually remain in the system for 
10 Myr in areas involved in the mean motion resonances (MMR) (e.g., 3:2 MMR ) 
with the inner companion. Moreover, we showed that the 2:1 and 3:1 resonances 
are on the fringe of stability, while the 5:2 MMR is unstable. Additionally, 
the 2:1 MMR marks out a remarkable boundary between chaotic and regular 
motions, inside, most of the orbits can survive, outside, they are mostly lost in the 
orbital evolution. In a dynamical sense, the most likely candidate for habitable 
environment is Earth-like planets with orbits in the ranges 0.8 AU < a < 1.0 
AU and 1.0 AU < a < 1.30 AU (except 5:2 MMR and several unstable cases) 
with relatively low eccentricities. The Trojan planets with low eccentricities and 
inclinations can secularly last at the triangular equilibrium points of the two 
massive planets. Hence, the 47 UMa planetary system may be a close analog to 
our solar system, bearing a similar dynamical structure. 

Subject headings: celestial mechanics-methods:n-body simulations-planetary systems 
stars:individual (47 UMa) 
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1. Introduction 



The main sequence star 47 UMa is of spectral type GO V with a mass of 1.O3M0. Butler 
& Marcy (1996) reported the discovery of the first planet in the 47 UMa system which has 
become one of the most eye-catching systems particularly after the subsequent release of an 
additional companion (Fischer et al. 2002). It is sometimes thought to be a close analog 
of our own solar system: for example, the mass ratio of the two giant companions in 47 
UMa is ~ 2.62 (see Table 1), as compared to that of Jupiter-Saturn (JS) of 3.34; and the 
ratios of two orbital periods are close to each other. Following the analogy, one may wonder 
whether there exists additional members in 47 UMa system. Multi-planet systems seem to 
be common in nature (Fischer et al. 2003; hereafter Paper I) more than ten of such systems 
have been detected in Doppler surveys to date. Moreover, if the solar system is a typical 
one, how are we to understand the analogous structures in other multi-planet systems, e.g., 
the presence of the inner low-mass terrestrial planets, the asteroidal belts or Kuiper belts 
(e.g. Spitzer Mission) around central stars? The pioneering numerical experiments (Adams 
& Laughlin 2003) showed that planet-planet scattering can yield planetary orbits with a ~ 1 
AU for a packed system consisting of 10 planets with widely spaced orbits. However, from an 
observational viewpoint, because of the variability of stellar atmospheres, the current ground- 
based observational precision ~ 2-3 m s^^ (Butler et al. 2003) indeed places a fundamental 
limit on the discovery of Earth- mass planets^ around 1 AU (even smaller), and the preference 
detection of planets with masses of several should benefit from future space projects that 
can carry out high resolution astrometric measurements (e.g., SIM and GAIA). Nevertheless, 
in the current stage, this leaves room for researchers to study the dynamical structure of 
such planets or their formation that may possibly reside in the systems, in advance of the 
observing missions. 



2. Formation scencirio 

The planetary formation scenarios may be replayed as follows: in the standard model 
of core- accretion scenario (Safronov 1969; Lissauer 1993), solid cores (or planetary embryos) 
rapidly grow larger by the accretion of kilometer-sized planetesimals via runaway growth 
(Kokubo & Ida 1998; Ida & Lin 2004), leading to the formation of the gaseous giant planets 
through the accretion of disk gas onto Earth-mass solid cores and more planetesimals into 



-'^however, the recent discovered Neptune-mass planets with close-in orbits (a < 0.1 AU) (Butler et al. 
2004; McArthur et al. 2004; Santos et al. 2004) are grouped as a new category linking giant gaseous planet 
and Earth-mass planet. 
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an envelope. When they formed outside the ice boundary ~ 3 AU (Ida & Lin 2004) in the 
system, two giant planets may migrate inward through disk-planet interactions (Ward 1997; 
Kley 2000; Nelson et al. 2000; Papaloizou 2003), and such migration process may not halt 
until the two planets are eventually captured into ~ 8:3 mean motion resonance (Laughlin, 
Chambers & Fischer 2002; hereafter Paper II), see also the 2:1 resonant capture scenario 
for two massive planets in GJ 876 (Lee & Peale 2002; Kley, Peitz, & Bryden 2004). In 
the migration scenario, the giant protoplanets may clear up wide gaps about their orbits 
(Goldreich, Lithwick, & Sari 2004) and cease accreting small bodies. But the giant planets 
may make the residual unaccreted planetesimals in the debris disk move inward or outward 
through gravitational scattering mechanism, even with the eccentricities and inclinations of 
these planetesimals excited by sweeping secular resonance (Nagasawa & Ida 2000), thus the 
orbits of the small bodies can imdergo mutual crossings and then they are either directly 
cleared up or collided into fragments in the post-formation stage. Another possible fate for 
these isolated members is that they may be captured into mean motion resonances with 
the giant planets under migration (e.g., the capture of the KBOs in the Neptune migration 
[Malhotra, Duncan, & Levison 2000]), and such resonant swarm may outline the structure of 
asteroidal belts in the system. In consequence, a dynamically quasi-stable configuration with 
spaced orbits may be achieved to support further secular evolution, which may determine 
the resulting destinies of those members, leading to the final assemblage of the planetary 
system (Lin 2004). 

Paper II and Gozdziewski (2002) studied the long-term stabihty of 47 UMa and pointed 
out that the secular apsidal resonance can help stabilize the two giant planets in an aligned 
configuration with the libration of their relative periapse longitudes (Ji ct al. 2003), then 
the eccentricities are well maintained to free from larger vibrations due to this mechanism, 
as a result, this system can even survive for billion years (Barnes & Quinn 2004). Several 
pioneer works were concentrated on the structure of the system and presented a preliminary 
understanding of this issue. Jones, Sleep & Chambers (2001) investigated the existence of 
Earth- mass planets in the presence of one known giant planet, and subsequently Paper II 
and Asghari et al. (2004) further studied the stability of massless test particles about the 
so-called Habitable Zones (HZ) according to some earlier solutions (Fischer et al. 2002), 
where the dynamical model was treated as a restricted multi-body problem. Nevertheless, 
as the terrestrial planets possess significant masses, they can interact with the two giant 
planets by mutual gravitation, which may result in secular effects for the planetary system. 
Accordingly, we should take into account the masses of terrestrial bodies in the model when 
exploring the dynamical architecture. In this paper, we performed extensive simulations 
to examine the dynamical architecture in both the HZ and extended areas, for Earth-hke 
planets (with masses from 0.1 to 10 M®) of 47 UMa with stable coplanar planetary 
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configuration, based on the best-fit orbital parameters given by Paper I. These new rehable 
orbital solutions are derived from additional follow-up observations, hence they can represent 
the actual motions of the system under study. On the other hand, as mentioned previously, 
the discovery of three close-in Neptune-mass planets demonstrates that it may be possible 
for less massive planets (~ M®) to move close to the star. Therefore, in the extended study, 
we also explored low-mass planets in the region 0.05 AU < a < 0.4 AU and we found that the 
secular resonance arising from the inner giant planet can render the eccentricity excitations 
for the Earth-hke planets (see §3.2). In addition, we also carried out two other runs to 
predict potential bodies in 47 UMa : In §3.1, to compare with the inner solar system, we 
further investigated the planetary configuration of 2 giant planets plus 4 terrestrial planets, 
and in §3.3 we examined the case of presence of Trojan planets with respect to the two giant 
planets. The results suggest that the 47 UMa system may have a similar architecture to the 
solar system. 

In the simulations, we use an N-body codes (Ji, Li & Liu 2002) of direct numerical 
simulations with the RKF7(8) and symplectic integrators (Wisdom & Holman 1991). We 
always take the stellar mass and the minimum planetary masses from Table 1. The adopted 
time stepsize is usually ~ l%-2.5% of the orbital period of the innermost planet, which 
is sufficiently small for the integration. Additionally, the numerical errors were effectively 
controlled over the integration timescale, and the total energy is generally conserved to 10~^ 
- 10~^ for the integrations. Our main results now follow. 

3. Simulations Results 

3.1. 2 Giants plus 4 terrestrial planets 

The announcement of the fourth planet in the 55 Cancri system (McArthur et al. 2004) 
suggests that multiple planetary systems resembling our solar system could be quite common 
in the galaxy. In the first runs, we examine the configuration consisting of 2 giant planets 
(2G) and 4 terrestrial planets (4T) to study the coexistence of multiple bodies. Specifically, 
this means that we directly place Mercury, Venus, Earth and Mars into the 47 UMa system 
to simulate "the inner solar system", where the orbital elements for above terrestrial planets 
are calculated from JPL planetary ephemerides DE405 at Epoch JD 2448750.9 corresponding 
to the outer companion (see Table 2), e.g., the semi- major axes are respectively, 0.387, 0.723, 
1.00 and 1.523 AU. The giant planets are always assumed to be coplanar in the simulations, 
thus the inclinations for 4T refer to the fundamental plane of the 2G's orbits. In Figure la, 
the orbital evolution of 4T is shown, where Q — a(l + e), q — a(l — e) are, respectively, 
the apoapsis and periapsis distances. The simulation results indicate that the orbits of 
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these low-mass planets become chaotic after 2 x 10^ yr: Mars leaves its initial orbit in less 
than 10'' yr and its Q can amount to ~ 10'^ AU at t > 1.5 x 10^ yr. Such dynamical 
instability stems from strong perturbation by 47 UMa b (see also §3.2.4), and the Hill radius 
is Rh = [Mi/(3Mc)]^''^ai (Mc, the stellar mass; Mi,ai, the planetary mass and semi-major 
axis of Companion B), then ^Rh — 0.6 AU, showing that the orbit of Mars is quite close to 
the outskirts of ?iRh sphere. In addition, as of its eccentricity grows, Venus begins to cross 
the orbits of Earth and Mercury at i > 2.0 x 10^ yr, which leads to eventual destruction of 
the system. 

However, as the two giant planets in 47 UMa are much closer to the central star than 
the JS-pair to the Sun, for perfect analogy with the inner solar system, the semi-major axes 
of 4T should be shifted to 0.15, 0.29, 0.40 and 0.61 AU. Bearing this in mind, we restarted 
a new run for the 2G-4T system, where we adopted the rescaled semi-major axes for 4T 
together with all the other initial values in Table 2. In this numerical experiment, we found 
that the 2G-4T system can be dynamically stable and can last at least for 20 Myr (sec Figure 
lb), where the time behaviors of Q and q for 4T show regular motions that their semi-major 
axes and eccentricities do not dramatically change in their secular orbital evolution. More 
terrestrial planets can be possibly created in the later stage of planetary formation. Chambers 
(2001) studied the accretion of ~ 150 planetary embryos with lunar-to-Mars masses in giant 
planets systems and recovered 4T similar to those in the solar system. However, Levison & 
Agnor (2003) underhned that the population and masses of the resulting terrestrial planets 
may be affected by the giant planets, because the growth and evolution of the planetary 
embryos are determined by the possibly experienced dynamical mechanisms (e.g., secular 
resonances) . 



3.2. Terrestrial planets in Habitable zones 

The Habitable Zones are generally conceived as places where the biological evolution of 
life is able to develop on planetary surfaces in environment of liquid-water, subtle temper- 
ature and atmosphere components of CO2, H2O and N2 (Kasting et al. 1993); at the same 
time, the planetary habitability is also related to the stellar luminosity and the age of the 
star-planet system (Cuntz et al. 2003). The HZ could be considered to be centered at ~ 1 
AU (Mc/Mq)^. For 47 UMa, the inner and outer boundaries of HZ range from 0.7 AU to 
1.3 AU (Mcnou & Tabachnik 2003), however, in our practical simulations, we extended the 
HZ to other areas for the purpose of a more comprehensive study. 

In the second series of runs, we extensively investigated the case of two giant companions 
with one terrestrial planet in the HZ. The mass of the assumed terrestrial planet ranges 
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from 0.1 to 10 M®. And the adopted initial orbital parameters are as follows: numerical 
scanning was carried out for [a, e] space by direct integrations, where the low-mass bodies 
were placed at equal intervals of 0.01 AU for 0.05 AU < a < 2.0 AU, the eccentricities 
were uniformly spaced every 0.01 for 0.0 < e < 0.2 (for 0.05 AU < a < 0.4 AU, where 
0.0 < e < 0.1), the inclinations are 0° < / < 5°, and the other angles were randomly 
distributed between 0° and 360°. Thus, over 3000 simulations were exhaustively performed 
for typical integration spans from 1 Myr to 10 Myr, totalling several 10^° yr. 

3.2.1. 0.05 AU<a<O.A AU 

In these runs, we explored the secular evolution of 385 "hot Earths" or "hot Neptunes" 
for a time span of 1 Myr. All the simulations arc dynamically stable for 10^ yr, and 96% of 
the orbits bear e final < 0.20. However, Figure 2a shows that the eccentricities for the bodies 
at ~ 0.30 AU are excited to ~ 0.40, where the secular resonance ui (4l".ll/yr) of the inner 
companion (similar to 1/5 for Jupiter) is responsible for the excitation of eccentricity. The 
debris disk at ~ 0.30 AU is also shown by Malhotra (2004), who presented similar results 
of the eccentricity excitation of massless bodies by nonlinear analytic theory for secular 
resonance. In addition, even more terrestrial planets with spaced orbits can simultaneously 
survive for longer lifetime (e.g., the second simulation of 2G-4T system in §3.1). Nevertheless, 
one can easily sec that this region is not a good location for habitability, due to extra high 
temperature. Recently, Narayan, Gumming, & Lin (2005) discussed the detectability of low- 
mass objects with orbital periods P < 10 days, and they estimated the present threshold 
for detecting planets with close-in orbits is > 10 — 2OM0, depending on the number of 
the observations and the precision of Doppler measurements. Indeed, the improvement 
of precision of the ground-based observations will lead to the discovery of additional low- 
mass planets (~ lOM®) in other known planetary systems (G. W. Marcy 2004, private 
communication) . 

However, it brings about great difficulty in catching the planetary formation for hot 
planets that move so close to their host stars. Although the gas giant planets may undergo 
inward orbital migration from several AU to the vicinity of host stars, Ida & Lin (2004) 
found that a large fraction (90%-95%) of the planets that have migrated to a < 0.05 AU 
must perish, because the tidal heating can make the planetary radius inflate (Gu, Lin, & 
Bodenheimer 2003) and the body may be directly engulfed by the star (Israclian ct al. 2001) 
or survive only as rocky "cores" with the gas envelope removed due to the mass loss of the 
planet through the material exchange at the Roche radius. Nevertheless, it is hopeful to 
detect such planets with a > 0.05 AU in future space missions (e.g., GOROT, KEPLER, 
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TPF), while the survival of the short-period terrestrial planets (Mardhng & Lin 2004) is also 
related to the relativistic potential of the star. 

3.2.2. 0.4 AU <a<l.Q AU 

We carried out 1260 integrations in this region for 5 Myr and we found that none of 
the orbits escaped during this time span and 94% of them were in the resulting e < 0.25, 
see the final status shown in Figure 2b. We find that the eccentricities of the orbits with 
0.70 AU < a < 0.78 AU can be pumped up and in the 2:9 MMR at ~ 0.76 AU, e can reach 
~ 0.90, indicating that there may exist a gap near this resonance. Most of the Earth-like 
planets about 1:4 MMR at ~ 0.82 AU move stably in bounded motions with low-eccentricity 
trajectories, except for two cases where the eccentricities eventually grow to high values. 
Paper II pointed out that the secular resonance V2 arising from the outer companion (similar 
to for Saturn) can remove the test bodies. Would the V2 also influence the Earth-hke 
planets in this system? Nevertheless, we did not find this mechanism at work at about 
~ 0.85 AU (see Paper II) when we examined the results, because the terrestrial planets 
under study that all bear finite masses that may change the strength of this resonance; on 
the other hand, the location of the secular resonance is changed due to the orbital variation 
of the outer companion. For a terrestrial planet with a mass of lOM®, the region for z/2 
secular resonance is now shifted to ~ 0.70 AU, where two eigenfrequencies for the terrestrial 
body and outer giant planet given by the Laplace-Lagrange secular theory are, respectively, 
21l".37/yr and 225".48/yr. This indicates that both planets almost have the same secular 
apsidal precession rates in their motion. At the new location, the V2 resonance, together 
with the mean motion resonance, can work at clearing up the planetesimals in the disk (see 
Fig. 2b) by the excitation of the eccentricity; qualitatively, our results are in accord with 
those of Paper II. 

While the inner edge of HZ marks out a narrow unstable area, it is very possible to 
discover Earth-like planets in this wider area, 0.4 AU < a < 1.0 AU, in future surveys, 
which are to be the best candidate of habitable places for biological evolution of intelligent 
beings. 

3.2.3. 1.0 AU <a<l.?> AU 

There were 630 simulations in this region for 10 Myr and we found that 88% survived 
the integration, confirming the results given in Paper II that most of the test particles with 
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a < 1.3 AU can eventually remain in the system. Here, for the 3:1 resonance at ~ 1.0 AU, 
Figure 3a shows that there are stable orbits in the zones about 1 AU with e < 0.1, which 
agrees with the work by Rivera & Haghighipour (2004) who showed that a test particle 
can last 100 Myr at 1 AU in 47 UMa; while for 0.1 < e < 0.2, the orbits tend to be in 
unstable state owing to the excitation of the eccentricities. The simulations may imply that 
the Earth- mass planets near 3:1 resonance are possibly on the edge of stability. However, the 
previous studies on this system showed that the 3:1 resonance is a gap with no survivors. Let 
us mention that such differences may arise from the adopted initial planetary configurations, 
and here we adopt the rehable best-fit orbital solutions given in Paper I that can describe the 
exact motions for the two giant planets. Hence, a comparative run was carried out to examine 
this, again we ran 630 simulations for 10 Myr but with the earlier orbital elements for the two 
massive planets to reproduce the previous results at the 3:1 resonance. Our results with the 
earlier data show that most of the Earth-mass planets about 3:1 resonance are unstable for 
the investigated time, and their eccentricities can be pumped up to ~ 1 through resonance; 
besides, the inclinations are excited to high values ranging 90° to 180°, indicating that the 
orbits of the Earth-mass planets become retrograde in the dynamical evolution and cross 
those of the prograde giant planets before they terminate their dynamical lifetimes. Thus, 
we may safely conclude that the stability of the terrestrial planets is dependent on the initial 
planetary configuration. 

In Fig. 3a, a narrow unstable stripe appears at the 5:2 MMR at ~ 1.13 AU, although 
several of them can be luckily left behind, most of the Earth-size planets are removed at ~ 1 
Myr due to the perturbation of 47 UMa b, and in this sense it is analogous to the situation 
in the solar system. However, a wider area between 3:1 and 5:2 MMR is assumed to be a 
qualified candidate habitable environment where the Earth-mass planet will not encounter 
the problem of dynamical stability, and this is almost true for the region (1.13 AU, 1.30 
AU) with e < 0.1, except several unstable islands near 7:3 MMR at 1.18 AU. The smaller 
eccentricity (near-circular orbits) may not cause dramatic variations of temperature on the 
planet's surface, so favoring habitabihty. Therefore, in a dynamical sense, if the 47 UMa 
system can be adopted as a candidate target for SIM, it is also possible to detect other 
Earths with stable orbits about 1 AU. 

3.24. 1.3 AU <a<1.6 AU 

We performed 651 simulations for 10 Myr, and found the dynamical structure in this 
regime to be quite complicated: 14% of them can finally survive for this time span, and 86% 
are lost by ejection into hyperbolic trajectories, indicating the chaotic nature for these bodies. 
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In Figure 3b, we can notice that the 2:1 MMR region is at ~ 1.31 AU and also close to the 
outer edge of HZ, and the orbits with 0.0 < e < 0.10 are unstable, while for 0.10 < e < 0.20, 
there are stable islands where fictitious planets can remain in bounded motions in the final 
system. We observe that the 2:1 resonance marks out a remarkable boundary between chaotic 
and regular orbits, indicating that orbits with a < 1.31 AU can have much larger surviving 
rates than those of a > 1.31 AU. However, there are wider stable region about 9:5 MMR 
at ~ 1.40 AU for low eccentricities 0.0 < e < 0.05. Most of the unstable orbits are in the 
region 1.43 AU < a < 1.56 AU, using resonance overlapping criterion (Murray & Dermott 
1999), the separation in semi-major axis Aa fa 1.3(Mi/Mc)^^^ai = 0.496 AU, thus, the inner 
boundary Rq for 47 UMa b is at ~ 1.58 AU, and the orbits in this zone become chaotic 
during the orbital evolution because the planets are both within ^Rh and also close to Rq- 
And the characterized ejecting time r ~ 1 Myr, which means the apparent gap (e.g., 5:3 
MMR at ~ 1.48 AU) in the inner belt, except for several stable islands. Another possible 
population for terrestrial planets is located at 3:2 resonance ~ 1.59 AU for 0.04 < e < 0.20, 
and 18 small bodies can last for 10 Myr and confirm the results of Paper II. The 3:2 MMR 
zone in 47 UMa is reminiscent of the Hilda asteroids in the solar system moving in a stable 
region. 



3.2.5. 1.6 AU < a <2.0 AU 

840 Earth-like planets are placed in this region for integration for 10 Myr, and the 
simulation revealed that 98% are removed from the system within a typical ejection time 
T < 3 X 10^ yr, which is much shorter than for 1.3 AU < a <1.6 AU, implying a thoroughly 
chaotic situation in this area. It is not difficult to understand that these terrestrial planets 
are entirely thrown out by the gravitational influence of the inner giant planet, using the 
SRh stability criterion. 



3.3. Trojan planets 

In the solar system, at present there are 1690 Jupiter Trojans^ (1062 preceding and 
628 trailing Jupiter) located at L4 and L5. Hence, a fascinating issue is that whether there 
exist "extrasolar Trojan planets" in other planetary systems, and if such planets can occur 
in the system harboring only one gas giant with a ~ 1 AU (Dvorak et al. 2004), it may 
essentially help understand their presence in the HZ. Laughlin & Chambers (2002) found 



^http://cfa- www.harvard.edu/iau/lists/JupiterTrojans. html 
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the equilateral configuration consisting of a star and two equal-mass planets can be linearly 
stable for mass ratios less than 0.03812 and showed that a pair of Saturn-mass planets can 
sustain the stability in the 1:1 orbital resonance. 

In the third series of runs, we let the experimental planets with masses of 0.1 to 
10 initially orbit about the triangular equilibrium points of the two giant companions, 
for each giant planet there are 100 orbits distributed near L4 (or L5) point. Then, we 
implemented 400 dynamical simulations to explore their secular evolution on the timescale 
of lO'^ yr. The starting conditions for the Trojan bodies arc: 0.00 < e < 0.10, 0° < / < 1° 
and |A — Ai,2 — 60°| < 1° (where A, and A1.2, respectively, the mean longitudes of the Trojan, 
the inner and outer companions). For the inner planet of 47 UMa b, the survival rates for 
the low- mass bodies about L4 and L5 points are, respectively, 96% and 98%, which is quite 
similar to the Trojans and Greeks for Jupiter; while for the outer planet (47 UMa c), the 
survival proportion goes down to 80% and 70% near L4 and L5 region, as a comparison 
Saturn has no co-orbital bodies. Obviously, the star-planet mass ratios n ~ 10~^ < 0.0385 
(Murray & Dermott 1999) both satisfy the linear stability condition for the triangular points. 
In addition, it is also suggested that the stable triangular region for the inner planet is wider 
than that of the outer planet, and this can be easily understood as this width depends on 
the value ~ //^/^ of the tadpole region. 

Figure 4 displays the long-term orbital evolution for the low-mass planets moving about 
L4 {black line) and L5 {yellow line) points for 47 UMa, where a and e both perform small 
modulations for 10^ yr, and A — Ai 2 also librate about 60° and 300°, respectively, with low 
amplitudes for the same time span. 

4. Summciry and discussions 

In this work, we have systematically studied the existence of Earth-like planets in the 
region for 0.05 AU < a < 2.0 AU for 47 UMa by numerical simulations. In addition, we also 
investigated the packed system and Trojan planets in this system. We now summarize the 
main results as follows: 

(1) The "hot Earths" can dynamically survive for 10® yr for 0.05 AU < a < 0.4 AU, and 
they are probably detected by transit time variations (Agol et al. 2005; Holman & Murray 
2005) via their interaction with the transiting planet'^. These techniques may become a 



^Laughlin G. and his collaborators are now maintaining the website (see http://www.transitsearch.org/) 
to search for transiting "hot Jupiter" planets through the worldwide cooperative observational cam- 
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vital and effective observational strategy to discover more transiting objects around main 
sequence stars. 

(2) The Trojan planets with low eccentricities and inclinations can last out at the trian- 
gular equilibrium points of the two massive planets of 47 UMa. In this sense, it is analogous 
to the case in solar system. Nevertheless, if they really exist, the formation of these bodies 
is still a mystery and needs further study. 

(3) The Earth-like planets can eventually remain in the system for 10 Myr in the areas 
associated with mean motion resonance (e.g., 3:2 MMR) with the inner companion. We 
also showed that the 2:1 and 3:1 resonances could be on the fringe of stabihty, but the 5:2 
MMR is unstable and the bodies can be ejected as "extrasolar comets" . And this may sketch 
out an asteroidal belt structure similar to the solar system. Moreover, the 2:1 MMR (near 
the outer boundary of HZ at 1.30 AU) marks out a significant barrier between chaotic and 
regular motions, implying that a large fraction of the orbits inside this resonance can be 
survive, while most of them are lost in the simulations outside the 2:1 resonance. Again, 
considering the inner boundary of HZ and dynamical stability, we point out that the most 
likely candidate for habitable environment is terrestrial planets with orbits in the ranges 0.8 
AU < a < 1.0 AU and 1.0 AU < a < 1.30 AU (except 5:2 MMR, and several unstable cases) 
with low eccentricities (e.g., 0.0 < e < 0.1). However, in our own solar system there are no 
terrestrial planets from the 1:4 MMR out to Jupiter, although there are stable orbits there. 
This may suggest that although some orbits arc stable, conditions arc such that terrestrial 
planets cannot form so close to giant planets. Perhaps this is because runaway growth is 
suppressed due to the increased eccentricities from the perturbations of the giant planet. 
In 47 UMa, the corresponding region runs from 0.82 AU on out (see Figure 2), almost 
completely covering the HZ. Hence, it would be reasonable to conclude that the only proper 
place to find habitable planets in this system would be at about 0.8 AU. But this should 
be carefully examined by forthcoming space measurements (e.g., SIM) capable of detecting 
low-mass planets. 

In a word, we can see that the 47 UMa planetary system may bear a similar dynamical 
structure to our solar system, in that it can also own several terrestrial members resembling 
the inner solar system. A comparative study can be performed in other planetary systems 
with two giant planets (Erdi et al. 2004; Raymond & Barnes 2005) to explore whether Earth- 
like planets can exist there, and to locate less massive planets (Malhotra 2005, in preparation) 
in general planetary systems. The dynamical habitability is relevant to both mean motion 



paign. Besides, the other active group from OGLE project have detected several transiting planets 
(http://sirius.astrouw.edu.pl/~ogle/) in the observations. 
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resonance and secular resonance for a given system. Moreover, the formation of the terrestrial 
planets in the HZ is difficult in the planetary systems that contain a giant planet with a 
moderate-size eccentric orbit (Veras & Armitage 2005), because most of the initial material 
within HZ was depleted by orbital crossings; in addition, the accretion simulations (Paper 11) 
show that the massive planetary embryos in the HZ should be created prior to the formation 
of giant gaseous planets. Hence, more efforts should also be directed to the understanding 
of the formation of habitable terrestrial planets. 
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Table 1. The orbital parameters of 47 UMa planetary system (adopted from Fischer et al. 

[2003]). The stellar mass is l.OSM©. 



Parameter 


Companion B 


Companion C 


Msmi{Mjup) 


2.86 


1.09 


Orbital period P(days) 


1079.2 


2845.0 


a(AU) 


2.077 


3.968 


Eccentricity e 


0.05 


0.00 


uj{deg) 


124.3 


170.89 


Periastron Time (JD) 


2452374.1 


2448750.9 



Table 2. The orbital elements for 4 terrestrial planets at JD 2448750.9 (Prom DE405). 



Planet 


a(AU) 


ecc 


I 


n 


CO 


M 


Mercury 


.387 


.2056 


7.00 


48.34 


29.10 


260.40 


Venus 


.723 


.0067 


3.39 


76.70 


55.06 


253.70 


Earth 


1.000 


.0163 


0.001 


255.71 


207.06 


123.84 


Mars 


1.523 


.0935 


1.850 


49.57 


286.48 


355.18 
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1000 : 

800 : 
600 : 
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Time(yr) 



Fig. 1. — The orbital evolution of 4 terrestrial planets. Left: chaotic case with a from Table 
2, Mars leaves its initial orbit at i < 10^ yr, and Venus begins to cross the orbits of Earth 
and Mercury at t > 2.0 x 10^ yr. Right: with rcscaled semi-major axes, 4T remain regular 
motions that a and e both perform slight vibrations for 20 Myr. 
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0.05 0,10 0.15 0.20 0.25 0.30 0.35 a(AU) 0.4 0,5 0,5 0,7 0,8 0,9 a(AU) 



Fig. 2. — The contour of status of the final eccentricities for Earth- hke planets, the vertical 
axis for the initial e. Left: 0.05 AU < a < 0.4 AU for 1 Myr. Notice z/i secular resonance at 
~ 0.30 AU pumps up the eccentricities. Right: 0.4 AU < a < 1.0 AU for 5 Myr. The e of 
the orbits with 0.70 AU < a < 0.78 AU can be excited and in the 2:9 MMR at ~ 0.76 AU, 
e can reach ~ 0.90. 
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1.00 1.05 1.10 1.15 1.20 1.25 a(AU) 1.30 1.35 1.40 1.45 a(AU)1.50 1.55 1.60 



Fig. 3. — The surviving time for Earth-like planets for the integration of 10 Myr, the vertical 
axis is the same as Fig.2. Left: for 1.0 AU < a < 1.3 AU, see the gap for the 5:2 MMR at 
~ 1.13 AU. Right: for 1.3 AU < a < 1.6 AU, a population of the terrestrial planets is about 
9:5 MMR at ~ 1.40 AU for low eccentricities (see texts for details). 
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Fig. 4. — The secular evolution for the low- mass planets moving about L4 {black line) and 
L5 {yellow line) points for 47 UMa, where a and e both perform small modulations for 10"^ 
yr, and A — Ai,2 also librate about 60° and 300° with low amplitudes, respectively. Left: for 
inner planet. Right.ioi outer planet. 



